Self-passivating edge reconstructions of graphene 
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Planar reconstruction patterns at the zigzag and armchair edges of graphene were investigated 
with density functional theory. It was unexpectedly found that the zigzag edge is metastable and a 
planar reconstruction spontaneously takes place at room temperature. The reconstruction changes 
electronic structure and self-passivates the edge with respect to adsorption of atomic hydrogen from 
molecular atmosphere. 
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Carbon is one of the most prominent elements in 
nature, vital for biology and life. Although macro- 
scopic carbon has been important since ancient times pQ, 
only modern materials design, utilizing nanotubes[2j [3], 
fullerenes[4] and single graphene sheets [5], fully attempts 
to use its flexible chemistry. In applications for nanoscale 
materials and devices, it is often the atomic and elec- 
tronic structure of boundaries and surfaces that is re- 
sponsible for mechanical, electronic and chemical prop- 
erties. 

Since the properties of nanomaterial depend on the 
precise atomic geometry, its knowledge is crucial for fo- 
cused preparation of experiments and for worthy theo- 
retical modeling. Only this enables the further develop- 
ment of nanoelectronic components, nanoelectromechan- 
ical devices and hydrogen storage materials [3j [6], or the 
usage of carbon in compound designs [7 . 

The importance of precise geometry is emphasized in 
low-dimensional systems. The strong correlations are 
known to bring up novel phenomena [8 , and such should 
be expected also for the quasi-one-dimensional edges of 
graphene. The edge chemistry plays even crucial role in 
the catalyzed growth of carbon nanotubes[9l [TU], Specif- 
ically, as two-dimensional carbon has the honeycomb lat- 
tice, edge behaviour ultimately boils down to the proper- 
ties of graphene edges. Hence it is relevant to explore dif- 
ferent graphene edge geometries and their chemical prop- 
erties beyond the standard zigzag and armchair ones. 

This relevance is evident from the abundant litera- 
ture. The electronic properties of graphene as well as 
carbon nanotube armchair and zigzag edges have been 
studied extensively [TTJ [12], often in connection to nan- 
otube growth [9l [13] or the so-called electronic "edge 
states" [14j [TSJ [16] . There has been experimental and the- 
oretical work done even on the reconstruction of graphene 
edges, but they have differed from the basic reconstruc- 
tion patterns studied in this work. They have involved 
either edge roughness [15] or more dramatic folding of the 
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FIG. 1: The geometries of graphene edges, (a) reconstructed 
zigzag (zz(57)), (b) armchair (ac), (c) reconstructed arm- 
chair (ac(677)), (d) zigzag (zz), and (e) pentagonal armchair 
(ac(56)) edge. The numbers in parentheses denote the num- 
ber of vertices in edge polygons. Some bond lengths (in A) 
and bond angles are shown on the right: the bond angles are 
a = 143°, (3 = 126°, 7 = 148° and 8 = 147°. All geometries 
are strictly in plane. 



edge into a loop [17]. 

The edges discussed in this work are shown in fig- 
ure [l] and were investigated by modeling an infinitely 
long carbon nanoribbon of a given width (see Methods). 
Tight-binding method [18 was used to explore number 
of other edge candidates, but density-functional analysis 
only for the relevant ones is reported here. The most 
important edge is zz(57), a reconstruction of zigzag edge 
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where two hexagons transform into a pentagon and a hep- 
tagon, like an edge cut through a Haeckelite structure or 
a line of Stone- Wales defects [T9]. The edge ac(677) is a 
reconstruction of the armchair edge where two separate 
"armrest" hexagons merge into adjacent heptagons by 
Stone- Wales mechanism. The pentagonal reconstruction 
ac(56) of armchair has a slightly different nature, since it 
requires the diffusion of carbon atoms from distant "arm- 
rests" to "seat" positions. 

Let us start analysis by looking at edge energy £ e dg e , 
which is calculated from the total energy of the graphene 
ribbon 
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where N is the number of carbon atoms, L the total 
length of edges (twice the length of simulation cell), and 
s gr = 7.9 eV is the cohesion energy of graphene. The 
edge energies converge rapidly as shown in figure [2] and 
justify the reference to (semi-infinite) graphene. The en- 
ergy of armchair edge is 0.33 eV/A lower than for zigzag 
edge, in accord with a similar value for nanotubes[T0]. 
However, the principal result is that by edge recon- 
struction zigzag may lower its energy by 0.35 eV/A. 
This implies that the reconstructed zigzag is the best 
edge for graphene. To the best of our knowledge, this 
is the first report on the metastability of the zigzag 
edge, which is surprising in view of the abundant litera- 
ture. Thermal stability of this novel reconstruction was 
also confirmed with tight-binding simulations [T8]. The 
ac(677) edge has only slightly higher energy than the 
armchair edge, whereas ac(56) reconstruction has the 
highest edge energy. We remark that the reconstruc- 
tions appear to be stable with respect to out-of-plane mo- 
tion, a situation somewhat different from small-diameter 
nanotubes[9l [T3]. 

These energetics can be understood by looking at the 
geometries of figure [I] Previous studies have shown 
that armchair has low energy due to triple bonds in 
the "armrests" [12] , as realized by comparing their short 
bonds (1.24 A) to the bond in acetylene (1.20 A). Zigzag 
does not have such triple bonds and ends up with strong 
and expensive dangling bonds. The reconstructed zz(57) 
has triple bonds (1.24 A) but also wider bond an- 
gles (143°) which reduces the hybridization energy cost. 
Triple bonds with wide angles can be observed also in 
ac(677) edge, but strain in other parts makes the recon- 
struction unfavoured. The ac(56) edge suffers from dan- 
gling bonds like zigzag, and additional high strain energy 
makes this reconstruction the most expensive one. Re- 
gardless, ac(56) edge has relevance during the growth of 
armchair edges [TO]. 

These observations are supported by hydrogen atom 
adsorption energies, shown in figure [3] The weak adsorp- 
tion for armchair (4.36 eV) compared to zigzag (5.36 eV) 
stems from the triple bonds in the armchair edge. Simi- 
larly the weak adsorption for zz(57) (3.82 eV) and ac(677) 
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FIG. 2: (color online) The edge energies of carbon nanorib- 
bons. Energies are plotted as a function of the ribbon width 
for the edges in figure [I] 
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FIG. 3: (color online) Hydrogen adsorption energies. The 
upper symbols correspond to e a ds with coverage of one hy- 
drogen per edge atom and lower faint symbols are shifted by 
subtracting the H2 binding energy; positive e a ds — Eh 2 means 
the hydrogen atom is more strongly bound to the edge than 
to H 2 molecule. 



(3.64 eV) witnesses the weakening of dangling bonds. 
The adsorption for ac(56) is large because of the strongest 
dangling bonds. For insight, the adsorption energies in 
figure [3] are replotted by subtracting hydrogen molecule 
binding energy Eh 2 — 4.58 eV. The resulting negative 
number for zz(57) means that the adsorption of hydrogen 
atom from H 2 molecule is not favored due to cost of H 2 
dissociation energy, unless the adsorption process should 
be complicated [20]. This amounts to the conclusion that 
the edge reconstruction chemically passivates zigzag edge. 
However, hydrogen adsorption for zigzag+hydrogen edge 
has yet smaller adsorption energy of 2.14 eV. 

The edge energetics are summarized in table [l[ Note 
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TABLE I: Summary of the edge and hydrogen adsorption 
energies. sl dge is the energy per edge atom, e a ds the hydrogen 
adsorption energy with full edge coverage, and s e dge+ads is the 
edge energy with hydrogen termination. Note that for ac(56) 
the edge atom density is A ac ( 56 ) = A ac /2. 
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zz(57) 


ac 


ac(677) 


zz 


ac(56) 


Sedge (eV/A) 


0.96 


0.98 


1.11 


1.31 


1.51 


zldge (eV/atom) 


2.36 


2.09 


2.30 


3.22 


6.43 


Sads (eV) 


3.82 


4.36 


3.64 


5.36 


5.58 


Sedge+ads (eV/A) 


0.34 


0.01 


0.45 


0.06 


0.74 



that the ordering of edges changes when expressed as 
energy per edge atom (e* edge = Sedge^ 1 ) due to differ- 
ent edge atom densities X ac — (2.13 A) -1 and X zz = 
(2.46 A) -1 . More interesting is to look at edges with 
hydrogen termination (Klein edge) [16 j. For this case the 
edge energy is e edge + ads = s edge -X-(e ads -E H2 /2), where 
the reference is to bulk graphene and H 2 molecules. The 
best edges in this case are normal armchair and zigzag 
edges because of dangling bonds. On the contrary, the 
weak dangling bonds and small adsorption energy causes 
high energies for zz(57) and ac(677) Klein edges. 

Let us now concentrate on the thermodynamic and 
electronic properties of zigzag edges. We would like to 
clarify an aspect which is ambiguous in the literature: 
there are two types of electronic zigzag edge states. The 
so-called "flat band" in figure [4^i comes from the bulk 
7r-electrons and is indeed localized at the edge. But the 
band due to dangling bonds is the one seen in the scan- 
ning tunneling microscope (STM) image and is located 
spatially even beyond the edge (figure^). In zz(57) for- 
mation of triple bonds is evidenced by the nearly isolated 
dimers in figure [4]i, and the reconstruction removes the 
dangling bond bands away from the Fermi-level by lift- 
ing the degeneracy almost by 5 eV. Hence for zz(57) the 
STM image shows only the "flat band" states. Because 
the dangling bond bands shift to elusive energies, also 
chemical reactivity reduces. 

In thermodynamic sense the spontaneous reconstruc- 
tion of zigzag into zz(57) should be possible, since the ac- 
tivation barrier from the zigzag side is only 0.6 eV, from 
the reconstructed side 2.4 eV. The G-mode vibration of 
graphene at 1580 cm -1 gives an attempt frequency of 
vq ~ 5 • 10 12 s _1 , and elementary approach yields the 
rapid rate vq • exp (— Es/ksT) ~ 4 • 10 2 s _1 at room 
temperature. 

The reconstructions predicted in this work are ex- 
pected to survive on graphite terraces due to the weak 
interaction (5.6 meV/atom) between the basal planes [2T]. 
Using appropriate sample preparation it should be thus 
possible to observe the reconstruction. STM images of- 
ten show irregular and blurred edges, yielding no atomic 
resolution, but at least for passivated edges armchair pre- 
dominance is claimed, in agreement with table [l[ 16] [22]. 



So far samples have been prepared intentionally with hy- 
drogen passivation during heat treatment [22} [23] . a sit- 
uation where reconstruction would not be favored. Al- 
ternate routes for observing the reconstruction would be 
the radial distribution function from neutron diffraction 
experiments without deuterium atmosphere [23 , or detec- 
tion of triple bond -spawned high-energy modes around 
2000 cm -1 with Raman spectroscopy. 

A further topic is the study of the conductance of var- 
ious graphene edges, particularly the zz(57) edge found 
here. The presence of the edge state around the Fermi 
level makes zz(57) ideal for conductance measurements, 
in contrast to armchair ribbons where the edge state is 
absent. This may render zz(57) as an interesting sta- 
ble model for quasi-one-dimensional carbyne with al- 




FIG. 4: (color) The electronic structure of zigzag and zz(57) 
edges, (a) and (b) shows the band structure for 34 A wide 
zigzag and zz(57) nanoribbons, respectively, with unit cell 
width of 4.9 A. Note that for zigzag this is twice the min- 
imum unit cell and the reciprocal space is thus only half of 
the normal representation. The dashed line is the Fermi-level. 
The bands coloured were identified directly by visual inspec- 
tion of the wave functions, (c) and (d) show the height pro- 
files of simulated scanning tunneling microscope images in 
constant current mode of the respective edges (height varia- 
tions > 2 A) , formed by integrating the electron density from 
occupied bands within 0.1 eV of the Fermi energy. The de- 
generacies at the gamma-point are 2 and 4 for the dangling 
bonds and the flat band, respectively. 
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ternating single and triple bonds. Furthermore, the 
novel thermodynamically and chemically stable recon- 
struction could play a role in formation of angular joints 
in nanoribbons [24] , closure of the ends of nanotubes after 
cutting [25] and any other system where graphene sheets 
are joined to produce systems with nano-scale morphol- 
ogy. 



Methods 

We used density- functional theory in conjunction 
with generalized gradient approximation for the ex- 
change correlation functional [26 and projector aug- 
mented waves [27] for the C(2s2p) electrons, as im- 
plemented in GPAW code using real- space grids [28 . 
Converged energies were obtained with grid spacing of 
0.2 A and 10 k-points in the periodic direction. In the 
perpendicular directions the system is not periodic and 
the space between the atoms and the wall of the sim- 
ulation cell was > 5.0 A. The energies were converged 
to ~ 10 -5 eV/atom and structures were optimized un- 
til forces were less than 0.05 eV/A. Our calculations 
agree well with previous relevant experimental as well 
as theoretical energetic and geometric properties [T2 ] |2T] . 
The activation barrier was calculated with nudged elastic 
band method [29, 30] by fixing the atoms beyond the first 
two zigzag-rows with unit cell of length 4.9 A. The con- 
stant current STM of Fig. [4] shows the height profile of 
electron density isosurface of the occupied states within 
~ 0.1 eV energy window below the Fermi- level. The iso- 
surface value corresponds to average density 2 A above 
the graphene plane (3 • 10 -5 electrons/A 3 ). 
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